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ABSTRACT 

We present high-resolution optical HST Space Telescope Imaging Spectrograph (STIS) observations 
made with two slits crossing four of the optically brightest starburst clumps in the vicinity of the 
nucleus of M82. These provide Ha kinematics, extinction, electron density and emission measures. 
From the radial velocity curves derived from both slits we confirm the presence of a stellar bar. We 
also find that the super star cluster M82-A1 has a position and radial velocity consistent with it being 
at the end of one of the unique X2 bar orbits formed by an inner Lindblad resonance. We derive a 
new model for the orientation of the bar and disc with respect to the main starburst clumps, and 
propose that clump A has formed within the bar region as a result of gas interactions between the bar 
orbits, whereas region C lies at the edge of the bar and regions D and E are located further out from 
the nucleus but heavily obscured. We derive extremely high interstellar densities of 500-900 cm~ 3 , 
corresponding to ISM pressures of P/k ~ 0.5-1.0 x 10 7 cm -3 K, and discuss the implications of the 
measured gas properties surrounding the nuclear star clusters on the production and evolution of the 
galactic wind. Despite varying pressures, the ionization parameter is uniform down to parsec-scales, 
and we discuss why this might be so. Where the signal-to-noise (S/N) of our spectra are high enough, 
we identify multiple emission-line components. Through detailed Gaussian line-fitting, we identify 
a ubiquitous broad (200-300 kms -1 ) underlying component to the bright Ha line, and discuss the 
physical mechanism(s) that could be responsible for such widths. We conclude that the evaporation 
and/or ablation of material from interstellar gas clouds caused by the impact of the high-energy 
photons and fast-flowing cluster winds produces a highly turbulent layer on the surface of the clouds 
from which the emission arises. 

Subject headings: galaxies: individual(M82) - galaxies: ISM - galaxies: starburst - galaxies: star 
clusters - galaxies: kinematics and dynamics - ISM: evolution 



1. INTRODUCTION 

Galaxies undergoing starburst events, triggered by in- 
teractions or mergers, are important objects to study in 
the local universe because of the insights they provide 
into the violent star formation processes that occurred in 
galaxies at much earlier epochs. M82 is the best-studied 
nearby (3.6 Mpc: [Freedman et al.lll99l) example of such 
a galaxy. The current starburst is thought to have been 
trig gered by a tidal intera ction with M81 a few x 10 s yrs 
ago (lYun. Ho. fc Lolll994l ), but investigations of this star- 
burst are hampe red by the galaxy's almost edge-on incli- 
natio n (i ~ 80°; iLvnds fc Sandagdll963t iMcKeith et~aTl 
|1995[) . Much of the starburst core therefore suffers from 
heavy extinction along the line-of-sight, so most observa- 
tions have concentrated on wavelengths beyond the vis- 
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ible domain (e.g. iRieke et al.lll980t lAchtermann fc Lacv 



1995| iSatvapal et al.lll995l . 119971 : IForster Schreiber et al 
20011 



O'Connell fc M angano' (|l978t ) identified a number of 



high surface brightness clumps or regions in M82, de- 
noted A, C, D and E that define the optical starburst 
core, and cover an area of ~500 pc in diameter. How- 
ever, the study of starbursts at resolutions attainable 
with the Hubble Space Telescope (HST) has revealed the 
existence of hyper-luminous compact star clusters, so- 
called super star clusters (SSCs), within sites of intense 
star- formation. With HST Planetary Camera imaging, 
lO'Connell et~aT1 (fl995h identified over one hundred can- 
didate SSCs_wjthin_the visible starburst, while more re- 
cently [MefoliOI] (|2005h catalogued a total of 197 young 
massive clusters in the starburst core with HST/WFPC2 
observations, and they associated 86 with region A alone. 
This incredible density of star clusters gives rise to a very 
unusual, highly energetic environment, and is very dif- 
ferent to anything we find in our own Galaxy or even in 

other less intense starburstS; 

In a companion paper ()Smith et all [2006. hereafter 
iPaper J ), we describe observations of an isolated SSC 
in region A, designated M82-A1, and find it to be sur- 
rounded by a surprisingly compact Hn region (r = 4.5 ± 
0.5 pc compared to the cluster radius, r = 3.5 ± 0.5 pc) 
at an unusually high pressure (P/k = 1-2 x 10 7 cm -3 K). 
The cluster H n region was found to have an ioniza- 
tion parameter of \ogU = —2. 24 ± 0.18, in excellent 
agreem ent with the IR study of IForster Schreiber et all 
(|2001h . These authors found, despite the seemingly 
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chaotic and rapidly varying nature of the ISM, a uni- 
form log U w —2.3 for scales ranging all the way from 
a few tens to 500 pc, leading them to suggest that a 
similar star-formation efficiency and evolutionary stage 
exists across the whole starburst region. Our observation 
of a similar ionization parameter extends the uniformity 
of the starburst conditions down to scales of only a few 
arsec s. An alternative model presented bv lDopita et al.1 
2002D associates a constant observed U with situations 
where dust effectively competes with gas in absorbing 
ionizing ph otons. W e briefly discuss this issue in Sec- 
tion [5J] In iPaperU we discussed how such a compact, 
high density Hll region should theoretically not exist 
at the estimated age of the cluster (6.4 ± 0.5 Myr) if 
it had evolved according to the standard model for a 
pressure-driven bubble. We found we could explain its 
existence through the high ambient pressures resulting 
in non-standard evolution. 

High ambient interstellar pressures have also been in- 
ferred from measurements of the size of supernova rem- 
nants (SNR) in the core of M82. Radio studies have 
shown that they are unusually compa ct (radii <4 pc; 
iMuxlow et all Il994t iPedlar et all Il999h . This suggests 
that, together with the results mentioned above, the den- 
sities and pressures are the dominant drivers of the H n 
region evolution and that pressure in the starburst core 
is a major factor influenc ing the starburst evolution (see 
also iRiebv fc Riekll200l . 

M82 hosts a stellar bar that is tho ught to have formed 
during its last encounter with M81 dTelesco et a.1.1 1 19911 : 
IWills et al.ll2000t iGreve et alj|2002t ). IWills et alT(|2000f ) 
developed a model of the M82 bar system by compar- 
ing the predictions of their analytical model to neutral 
(Hi, CO) and ionized ([Ne ii]12.8/im) gas observations. 
They found that the ionized gas is dominated by a faster 
rotating component that they identify with so-called X2 
bar orbit fa mily (thus corresponding to the ionized ring 
identified by lAchtermann fc Lacv|[l995D . whilst the neu- 
tral and molecular gas appears more consistent with the 
slower rotating cusped X\ orbit family (a^- orbits are asso- 
ciated w ith an inner Lindblad re s onanc e; lAthanassoulal 
Il992al) . lAthanassoula fc Bureaul (|1999f ) perform exten- 
sive hydrodynamical simulations of bars, concentrating 
specifically on the gaseous component, and find the gas 
in fact follows stream-lines that correspond only qualita- 
tively to the x\- and a^-orbits. Detailed inspection shows 
that there are differences arising from the fact that in re- 
ality the gas strea m-lines only l o osely resemble the stellar 
orbit families. In IGreve etail (|2002f) . K. Wills and her 
collaborators used the near-IR CaiiA854 2 stellar photo- 
spheri c absorption line observations from lMcKeith et al.l 
(|1993h to compare their bar model directly to stellar 
velocities rather than the previously used gas mea sure- 
ments. By making this comparison, IGreve et a.1.1 were 
able to estimate that the mass of the stars following x^- 
orbits is ^15 per cent of the total mass of the bar. 

The densely packed star clusters, high densities and the 
presence of a bar give rise t o unus ual, highly energetic 
ISM conditions. lLord et all (|1996[ ) developed a model 
for the ISM in the central 700 pc of M82 based on far- 
IR spectroscopic observations of fo rbidden-lines (see also 
iForster Schreiber et al1l2001l 120031) . They proposed that 
the observed emission can be explained as originating in 
multiple Hll regions and PDRs (photo-dissociation re- 



gions) mixed uniformly with the ionizing stars and clus- 
ters. They found that the clouds have a characteristic 
size, r c \ ~ 0.4-1 pc (with a molecular core and extended 
ionized envelope), density, n c ~ 10 4 cm -3 , and mass 
range, M c \ ~ 200-3000 M Q , and are highly pressured 
(P/k ~ 3 x 10 6 cm~ 3 K). These states are maintained 
by a hot (~10 6 K), diffuse surrounding medium of equal 
pressure, supported by the influence of a large number of 
SN shoc ks. In order to maintain the observed level of ion- 
ization, lLord et al.l found the average separation of the 
clouds and ionizing sources had to be of order their size 
(1-7 pc). This agrees well with the aver age cluster sepa- 
rations catalogued bv lMelo et al.l (12005 ) , and with high- 
resolu tion CO observations (Mao et al. 2000; Wci fi et all 



1 2 If) that also show that the clouds are partly disrupted 
and dissocia ted. By mo delling the observed forbidden- 
line profiles. lLord et all argued that the clouds are clus- 
tered in two 'hot-spots' of size 125 x 125 pc located at 
~±15" either side of the nucleus, representing a cross- 
section through the molecular and ionized gas torus 
found to surround the nuclear bar (lLarkin et al.l 119941: 
Shen fc Lolfl995l : lAchtermann fc Lacvl 119951: iWeifi et al.1 
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M82 exhibits one of the largest optically vis- 
ible outflows or 'superwinds' in the local Uni- 
vers e. The outflow is centre d on r egions A and 
C tiShopbell fc Bland-Hawthorn! [l998t lOhvama et all 
l2002f) . and is composed of a complex morphology of 
loops and filaments. Recent ground- and space-based 
narrow-band optical imaging have shown that, con- 
trary to previous thought, the M82 outflow is com- 
prised of many channels which seem to point ba ck to 
the individual starburst clumps (| Wills et al.l 119991 . Gal- 
lagher et al., in prep.), suggesting that the wind en- 
ergy is injected by multiple discrete sources (i.e. SSCs), 
rather than a monolithic starburst regio n (see e.g. 
iTenorio-Tagle. Silich. fc Munoz-Tun6nll2003f ). The study 
of the ISM environment in which SSCs exist (the star- 
burst core) and from which the superwind is driven, can 
hope to give insights into the initial conditions for galac- 
tic wind formation. 

Observations of the ionized gas in the central region of 
M82 at visible wavelengths with HST offer the potential 
of probing the starburst region on sub-arcsecond scales 
(O'.'l 1.8 pc) and further quantifying all of the above 
issues. In this paper, we present spectroscopy obtained 
with the Space Telescope Imaging Spectrograph (STIS) 
of the ionized gas near the nucleus, including regions A 
and C. We measure the kinematics, densities and ex- 
tinction of the ionized gas on parsec-scales in order to 
characterise the nature of the ionized gas in the M82 
starburst. 

The remainder of this paper is organised as follows: 
in Section [3J we describe the data and the reduction 
methods employed, including a description of how we ex- 
tracted spectra at regular intervals along both slits and 
fitted the resulting emission-line profiles. The results of 
our line fitting are described in Section^ together with a 
discussion of the kinematics and line dynamics. In Sec- 
tion [4] we present our analysis of the properties of the 
ionized gas, including our flux, extinction and electron 
density results; in Section [5l we go on to describe the 
physical state of the gas and compare our results to pre- 
vious models of the ISM. In Section [6] we analyse the 
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Fig. 1.— HST/ACS WFC F658N image with the prominent clumps A, C, D and E identified, together with the location of the SSC 
M82-A1 and the 2.2 )im "nucleus" (white cross). The location of the two 52 X 0.1 arcsec STIS slits are also marked, where portions of the 
2D slit images corresponding to the white shaded parts of the slits are shown in Fig. [2] The dashed line perpendicular to the slits marks 
the zero reference point used in the following plots. 



structure of the starburst and derive a new model for 
the orientation of the bar and disc with respect to the 
starburst clumps. In this section we also discuss what im- 
plications our results have for our understanding of the 
supcrwind. Our findings are summarised in Section [7J 

2. OBSERVATIONS AND ANALYSIS 

As part of a larger programme devoted to spectroscopy 
of SSCs in M82 (GO 9117; P.I. O'Connell), we obtained 
HST /STIS spectra for two slit positions crossing the 
galaxy's starburst core. Fig. [T] shows the position of the 
two STIS slits overlaid on an HST/ACS WFC F658N im- 
age (GO 9788; P.I. Ho). The slits (52 x 0.1 arcsecs) pass 
close to the "nucleus" (2.2 /im peak) and just above and 
below the dust lane that perpendicularly bisects the disc 
of M82. The first slit, which we will refer to as slit Al, 
intersects the core of region A and the southern edge of 
region C, and was centred on the brightest isolated clus- 
ter in region A (which we have designated M82-A1) with 
position angle PA = 229°. The secon d slit, sli t B2, was 
positioned on a bright cluster (B2-1; IPaper J ) towards 
the north-east of the starburst core with a position an- 
gle PA = 235°, and passes through portions of regions D 
and E. The coordina tes and exposure times for each slit 
position are given in iPaperll . 

2.1. Description of 2D spectra 

Spectra were taken with two gratings giving coverage 
from 2900-5700 A with the G430L and 6295-6865 A with 
the G750M setting with a spatial resolution of 0.05 arc- 
sec pixel -1 . A detail ed descri ption of the data reduction 
procedure is given in IPaper 1 The spectral resolution of 



the data was measured for the G430L data by taking the 
average FWHM of Gaussian fits to the unresolved H/3 
line at various points along the length of the slit where 
extended H/3 emission is observed. We found the resolu- 
tion to be 2.5 pixels (in good agreement with the quoted 
value of 2-3 pixels given in the STIS Data Handbook for 
extended objects), corresponding to 6.7 A in the G430L 
grating and 1.4 A in the G750M grating. Since variations 
in slit illumination and of the source structure at each 
point could affect the resolution slightly, it is possible 
that some of the narrowest lines we measure (uncorrected 
FWHM < 75 kms -1 ) may be unresolved. The material 
covered by the majority of the north-easternmost portion 
of both slits is very faint (see Fig. [TJ, and not detected in 
our obse rvations. The extinct ion i n this region as mea- 
sured by lSatvapal et all (|1995f ) and lAlonso-Herrero et al.l 
(|2003h may be very high, Ay ~ 10-15 mag. In spec- 
tra where we do detect emission, we observe the nebular 
emission lines of H/3, [Om]A5007, [N n] AA6548, 6583, Ha 
and [Sn]AA6716,6731. 

Fig. [2] shows the Ha and [N n] spectral region of the 
G750M 2D spectral images for slits Al and B2. Both 
2D slit images are able to share the same x-axis since 
their position angles (PAs) are sufficiently similar (scale 
is given in arcseconds offset from the position of M82-A1 
- marked with a dashed line on Fig. [T|) . The spatial ex- 
tent of the slit images is also restricted to show only the 
nuclear regions where emission is actually detected (rep- 
resented by the white shaded section of the slits in Fig.[T]). 
We see emission arising from both discrete sources and 
holes in the foreground dust screen. Bright emission from 
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Fig. 2. — STIS G750M two-dimensional spectral images showing the Ha and [Nnl emission lines over the portion of the slits covering 
the nuclear regions (equivalent to the white shaded parts of the slits shown in Fig. [lj. The lower plot shows the image from the Al slit, 
and the upper plot shows the B2 slit. In both cases, the x-scale is in arcseconds offset from cluster Al (for reference, this point is marked 
with a dashed line on Fig. [TJ and the brightness scale is in units of 10 — 16 erg cm -2 s -1 A -1 arcsec - 2 . The emission lines are identified, 
and the slit orientation is marked. Spectra were extracted using 10 and 20 pixel apertures along the length of both slits (see text); arrows 
marking particular spectra show the position of the Ho profiles shown in Fig. [4] (where the numbers correspond to the plot titles). 



the ambient interstellar gas is also present. 

Fig. [2] clearly shows how the emission line profiles vary 
continuously along the length of the slit. It would be 
ideal to track these variations on a pixcl-to-pixel basis, 
but in order to make accurate measurements of the gas 
conditions, the spectra need to have a sufficiently high 
signal-to-noise (S/N) ratio. Therefore we summed the 
spectra using multiple-pixel apertures. To measure the 
kinematics of the gas, observations of a strong emission 
line (e.g. Ha) are all that is required. However, to derive 
gas properties that require information from additional, 
fainter lines, a higher S/N was needed. To satisfy these 
two criteria, we extracted two sets of spectra from each 
of the 2D images. A high S/N set for each slit was made 
by stepping along in 20 pixel intervals (= l'.'O) , extract- 
ing apertures of width equal to this interval (hereafter 
referred to as set A20), and resulted in spectra from slit 



Al with an average S/N of 10 in H/3 (G430L grating) and 
20 in Ha (G750M grating). In order to maximise spatial- 
resolution for the kinematics measurements whilst main- 
taining sufficient S/N for accurate line fitting, we de- 
creased the extraction widths to 10 pixels (= Of! 5) sepa- 
rated by intervals of 10 pixels to create the second data- 
set (hereafter referred to as set A10). Extracted spectra 
were given sequential IDs starting from the north-east 
of each slit to identify them in later analysis. No back- 
ground subtraction was attempted for either of the slit 
Al or B2 datasets due to the bright and variable nature 
of the diffuse emission in the starburst core. 

Fig. [3] shows four representative spectra from regions 
A, C, D and E illustrating the variation between these 
clumps. The arrows on Fig.[2]mark the position that each 
spectrum shown in Fig. [3] was extracted from, where the 
numbers indicate the spectrum IDs referred to above, 



The environment in the starburst core of M82 5 



)2 - 86 



CO 

u 

CD 




)2 - 100 



150 - 



100 



6500 6600 6700 6800 
Al - 61 



6500 6600 6700 6800 

A1 - 86 



r\ 150 - 



100 - 




6500 6600 6700 
Wavelength (A) 



6800 



6500 6600 6700 
Wavelength (A) 



6800 



Fig. 3. — Representative STIS G750M spectra of regions A, C, D and E. The numbers in the plot titles correspond to the arrows in 
Fig. [2] The emission lines are identified on the lower-left plot only. 



and correspond to the individual plot titles in Fig. [3l 

2.2. Cluster distributions and approximate ages within 
the clumps 

The clumps of emission seen in Figs [T] and [5] represent 
a densely packed mixture of individual clusters and their 
surrounding Hn regions, ionized gas clumps and bright 
diffuse ISM seen through holes in the foreground dust 
screen, all with sizes and separations on the order of a few 
to a few tens of pc. The extent to which these phenomena 
dominate the physical nature of the spectra along the 
length of the slits varies dramatically. 

In the cores of the starburst clumps, the slit s cross 
a number of individual clusters (see e.g. iPaper J . figure 
1 inset). We have experimented with extr acting spec- 
tra clo se to individual sources catalogued by lMelo et al.1 
(|2005| ). but unfortunately the slits to do not pass directly 
through the centre of any isolated clusters, making the 
extraction of unique spectra difficult. Those that we did 
extract were of too low S /N to allow accurate model fit- 
ting. Measuring the equivalent width (EW) of the Ha 
line in these core regions can, however, give an approxi- 
mate indication of age trends in the stellar population by 
comparing to evolutionary synthesis models. It must be 
borne in mind, however, that since there is no way that 
we can measure the light from individual clusters, effects 
such as crowding, slit losses, and the presence of bright 
ambient diffuse emission will strongly affect the meaning 
of our EW measurements. 

By measuring the Ha EW along the central core re- 
gions of both slits and using the EW pred ictions of 
a STARBURST99 fSB99; iLeitherer et al.|[l999f) . 10 6 M 



model with a Kroupa IMF formulation, we find that the 
cluster ages in regions A, C, D and E lie between 5.5 and 
6.5 Myr (log(EW) « 1.5-2.5). If, however, a significant 
fraction of the Lyman continuum photons are absorbed 
by dust (jDopita et al.ll2002l l2006al . and discussed in Sec- 
tion l5.1[) . then cluster ages determined in this way will be 
overestimated. Nevertheless, these values are consistent 
with previous age estimates of clusters in the starbu rst 
core (jForster Schreiber et al.ll2003l ; ISmith et al.ll2006h . 

2.3. Line Profile Fitting 

The S/N and spectral resolution of the 10-pixel wide 
spectra are high enough to resolve multiple components 
in Ha in many regions along the slit. In general, a bright, 
narrow component is superposed on an underlying broad, 
fainter component, but obviously the identification of the 
broad component is dependent on the S/N of the partic- 
ular spectrum. The component profiles are best approx- 
imated by a Gaussian function since the line broadening 
mechanism is dominated by Doppler effects caused by 
the gas temperature or turbulent state. 

Examples of Ha line profiles from different parts of re- 
gions A, C, D and E are shown in Fig.[l]together with the 
Gaussian profiles needed to model the integrated shapes 
(see below) . Again, the position along the slit that each 
Ha profile has been extracted from is marked by num- 
bered arrows on Fig. [2] These were chosen to represent 
a range in the variety of line shapes seen: bright, nar- 
row emission lines with underlying broad components; 
double-peaked lines with no detectable underlying broad 
emission; and mixtures of the two. Where a broad com- 
ponent is present, it is often, but not always, rcdshiftcd 
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Fig. 4. — Example Ha line profiles for slit Af (upper six) and slit B2 (lower two): arrows on Fig. [2] indicate where the spectra are 
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compared to the narrower component. Double-peaked 
lines with velocity differences of ~50-100 kms -1 are seen 
in both the cores of region A and C. 

We have used the general-purpo se IDL-based curve- 
fitting utility, pan (Peak ANalvsis: lDimeoll2005| ). to au- 
tomate the Gaussian fitting process of each of the ex- 
tracted spectral lines. Briefly, pan embodies an inter- 
active environment for specification of a profile's initial 
guess parameters, and uses a \ 2 minimisation algorithm 
to optimise the fit to the data. A more detailed descrip- 
tion of the program and how we have customised it for 
our use is given in Westmoquette et al. (2007a). 

Each line profile was fit with the appropriate number 
of Gaussian components, giving a measurement of the ra- 
dial velocity, flux, and FWHM for each component. The 
lowest S/N lines could only be fit with a single Gaus- 
sian. However along most of the slit where emission is 
detected, the line profile shape distinctly contains more 
than one component. In cases where a double-component 
fit was most appropriate (determined from a comparison 
of the x 2 n t value to the single component result, and 
by visual inspection), the initial guess was always made 
setting the first Gaussian as the narrow component (here- 
after referred to as cl ) , and the second Gaussian as the 
broader component (hereafter c2); in the case of a double 
peak, the bluest component was assigned as cl and the 
redder one as c2. For a few cases a triple-component fit 
was necessary, and the additional component (c3) was al- 
ways specified with a initial guess assigning it to a supple- 
mental narrow line at the same wavelength as the main 
narrow line. This consistent approach helped limit the 
confusion that might arise during analysis as to which 
Gaussian fit belonged to which component of the line, as 
well as aiding the \ 2 minimisation process employed by 
pan. This assignment convention means that wherever 
an underlying broad component exists, we always refer 
to it as c2. 

3. KINEMATICS OF THE IONIZED GAS 

There have been a number of previous studies aimed 
at measuring the characteristics of the ga seous com- 
pone n t in the disc and wind of M82 (e.g . iGotz et al.l 
19901 iHeckman et all fl99(l iMcKeith "etall 1199.1 I1995T 
Shopbell fc Bland-Hawthornlll998f ). From these studies, 
it is apparent that the wind outflow, as defined by two 
well-separated (~300 kms -1 ) velocity components, is ob- 
served out to distances greater than 1 kpc from the nu- 
cleus along the mi nor axis. The disc of M8 2 is inclined at 
an ang le of ~80° (|Lvnds fc Sandagdll963t) such that the 
near side of the disc is projected on to the north-west 
side of the nucleus and the nuclear regions, including 
clumps A, C, D and E, are observed from underneath 
through the south ern side (although this has been dis- 
puted in the past; IWills et alJl2000f >. We therefore di- 
rectly view the roots of the wind near regions A and C. 
The t wo STIS slits are inside the 300 p c wind injection 
zone (IShopbell fc Bland-Hawthornlll998T ) where disc ma- 
terial is being entrained into the flow, and likely sample 
much of the inter-cluster material within the main star- 
burst clumps. 

3.1. Radial velocity variations 

The measured radial velocities of the Ha emission-line 
components from Gaussian fits to the A10 spectrum- 



sets are shown in Fig. [5j and are plotted with respect 
to the systemic velocity of M82 (f svs = +200 k ms -1 ; 
Ide Vaucouleurs et all fl99lt IMcKeith et all 119931 ). We 
have again plotted the datasets from both slits together. 
As described in the above section, in general the Ha line 
profile in the brightest parts of all four regions is com- 
posed of a bright, narrow component superimposed on a 
broad, fainter component with a FWHM ~ 100 kms -1 . 
At a few positions, we see a third component (c5). Error 
bars are shown for points to the north-east of M82-A1 
and for all c3 detections. For most of the other points, 
the error bars are smaller than the size of the plotting 
symbol. 

The brightest component, cl, is expected to trace the 
velocity of the densest (i.e. bulk of the) ionized gas in 
M82. From Fig. El the slit Al velocity curve for cl 
reaches a maximum velocity of «+100 kms -1 at an 
offset of ^+2", flattens out between —2" to —1" co- 
incident with the conspicuous dust-lane bisect i ng th e 
galaxy disc (see Fig. |T| also lO'Connell et alJ Il995f k 
then begins to rise again beyond —7". The velocity 
of M82-A1 (jPaper If ), together with two points either 
side of it, are redshifted by ^40 kms -1 from the gen- 
eral trend. On the opposite side, the rotation curve 
shown by cl reaches a maximum blueshifted velocity 
of 130 kms -1 then begins to decrease again beyond 
+20". At these offsets, however, the A2 slit is far- 
thest from the major-axis and begins to sample the in- 
ner regions of the wind - we return to this point be- 
low. We measure a cl Ha velocity gradient for the 
region nearest to the nucleus of 12 kms -1 arcsec -1 , 
in excellent agreem ent with the measurements of 
IMcKeith et~aTl (|1993h and IShopbell fc Bland-Hawthorn! 
( 1991 of 11 kms" 1 arcsec . The radial velocities of the 
cl gas in slit B2 follow the same velocity trend as the gas 
sampled along slit Al, but with a shallower gradient of 
7 kms -1 arcsec -1 . We discuss what these measurements 
tell us about the rotation curve in Section [5J 

The broad underlying component, c2, largely follows 
the cl velocity curve although there arc some devia- 
tions. In the core of region C, c2 is redshifted by up to 
60 km s -1 compared to cl. Here we detect a third compo- 
nent (c5), which is redshifted by a further 30-40 kms -1 , 
and may represent the front side of an expanding struc- 
ture and possibly the roots of the wind flow. Further 
along the slit, to the south-west of region C, c2 becomes 
blueshifted by ~40 kms -1 with respect to cl, and there 
are two detections of a 40-80 kms -1 blueshifted third 
component. Again, these could be signatures of expand- 
ing structures and/or ordered flows. An unusual situa- 
tion occurs at ~+4" offset (right in the core of region A), 
where cl is redshifted compared to its expected velocity, 
and it is c2 that instead conforms to the trend. Here we 
also detect a third component, blueshifted by ~80 kms -1 
compared to cl, which might explain the unexpected ve- 
locities if they are both two halves of an expanding shell 
where the faintest component is the nearest. In slit B2, 
we can only detect c2 in the brightest parts of regions D 
and E, and find that their velocities are similar in all cases 
to that of cl, except in the core of region E where there 
is evidence of a 20-40 kms -1 offset. Critically, however, 
none of the radial velocity differences measured equal the 
galactic wind velocity shifts of ~300 kms -1 . This evi- 
dence suggests that we are observing inside the injection 
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Fig. 5. — Heliocentric Ha radial velocities measured from the A10 spectrum sets, in units of kms 1 with respect to the systemic velocity 
of M82 (v Bya = +200 kms" 1 ); error bars are not shown because they are approximately the same size as the symbols used. The x-axis 
scale corresponding to both slits is measured in arcscconds offset from the position of M82-A1. The projected distance along the slit from 
the location of cluster Al to the nucleus is ~5". Measurements from both slits are shown on the same plot for ease of comparison; the 
approximate extent of regions A and C that coincide with slit Al, and regions D and E that coincide with slit B2 are marked, cl refers to 
the bright, narrow component, c2 to the broad component, and c3 to the longer wavelength component when the line is resolved into two 
narrow components. 



zone where organized flow velocities ar e low (consistent 
with lShopbell fc Bland-Hawthorn|[l998l) and/or that the 
flows at this point are oriented in the plane of the sky 
(i.e. in a poleward direction). 

3.2. Line widths and components 

Fig. [5] shows the FWHM of the individual Ha line 
components for both slits, corrected for the instrumen- 
tal resolution (1.4 A) but not thermal broadening. As 
mentioned in Section [21 small variations in the spectral 
resolution due to changes in the slit illumination/source 
structure could mean that the narrowest lines are only 
marginally resolved. The large scatter in cl FWHM val- 
ues (left-hand panel) across the two slits can be partly 
explained by a contribution from multiple unresolved 
kinematic components. We have attempted to show 
this by highlighting cl profiles from both slits that only 
required a single Gaussian to fit with circled symbols: 
the FWHMs of circled points are, in almost every case, 
broader than the narrow component (cl) of a multiple- 
Gaussian fit. This is particularly obvious in the —9" to 
0" region of slit Al where a second component cannot 
be fit at all due to the low S/N of the extracted spectra. 
It is very unlikely that these differences in width reflect 
real changes, so we conclude that multiple components 
are present everywhere, but the S/N is not sufficient in 
all cases to permit a more accurate fit to be made. This 
argument is strengthened by the observation that in ev- 
ery spectrum extracted from either slit where the Ha 



line has a high enough S/N, we see evidence of a c2- or 
c5-type component. 

To measure the following averages, we ignore all single- 
Gaussian-only fits, as these are assumed to be biassed by 
unresolved components. For slit Al, we find an average 
FWHM of ~60 ± 40 kms" 1 for cl (whereas the average 
width of single-Gaussian-only fits is — 160 ± 40 kms -1 ). 
For c2, we find an average FWHM of —210 ± 60 kms" 1 , 
but in a few cases the width of c2 approaches 300 km s" 1 . 
Where detected, the mean FWHM of c3 is —50 kms" 1 , 
which is consistent with the mean line width of cl. In 
general, the FWHM of the gas in slit B2 is narrower: in 
region D the width of cl is as low as 20-30 kms" 1 , with 
a corresponding c2 FWHM of 150-200 kms" 1 . In parts 
of region D and the area between regions D and E, only 
one component can be fit (shown by circled symbols), 
with a width of —100-140 kms" 1 . The line widths in 
region E are similar to region D, and in both regions c2 
is never >200 kms" 1 . 

3.2.1. Line broadening mechanisms 

Fig. [6] shows that all but a few line components are fit 
with Gaussians with widths > 30 kms" 1 . Interestingly, 
the widths are systematically narrower in regions D and 
E (even including those that could only be fit with a 
single component), possibly indicating a more quiescent 
gas state. We will now explore the possible explanations 
for the observed widths of the narrow (cl) and broad 
(c2) lines in turn. 
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Fig. 6. — Ha FWHMs measured from the A10 spectrum sets for cl (left-hand panel) and c2 &; c3 (right-hand panel). The units are in 
kms -1 corrected for instrumental broadening, and the x-axis scales are measured in arcseconds offset from the position of M82-A1. Circled 
points are those where no more than one component was required to fit to the line. 

Supersonic line widths have been observed in many 
intense st ar-formation sites both in nearby galaxies 
(30 Dor: IChu fc Kennicuttl Il994t iMelnick et all 119991 ; 



NGC 604: |Yang et alJ Il996t NG C 2363: kov etall 



[l^lGonzaiez-Delgado et aljil994D and in more distant 
dwarf galaxies (llzotov et al.lll996l iHomeier fc Gallagher! 
19991 IMarlowe et all 119951: iMendez fc Estebanl 119971; 



Sidoli et all 120061: IRozas et al.l l2006fh In these stud- 
ies, gravitational br oadening through virial motions of 
ioniz ed gas clouds (jMelnickl Il977t iTerlevich fc Melnickl 
1 19811) . and multiple unre solved kinematical comp onents 
along the line-of-sight (|Chu fc Kennicuttl 11994) have 
been proposed as explanations for the broadening of the 
brightest component (what we refer to as cl ) . A number 
of other m echanisms propo sed, such as dust or electron 
scattering (jRov et all 11993) . can be immediately disre- 
garded on physical grounds. 

The fact that we rarely measure FWHMs < 30 kms -1 
across both slits indicates a global origin for a 
large fraction of the observed widths. A clear 
lower-limit to the width of cl observed in some 
studies of similar environments (iMuhoz-Tuhon et al 
1996t|Martmez-Delgado et al.l[2007HWestmoquette et al 

im 



2007allbf ) supports this explanation. The only two broad- 



ening mechanisms that can work on a global scale are 
gravitationally induced virial motions and the stirring ef- 
fects of wide-scale, intense star-formation on the ambient 
ISM. We therefore propose that the turbulent motions 
induced by these two mechanisms provide the observed 
base level of line-broadening present over the whole of 
the region observed. However, neither of these globally 
determined values can account for the large local varia- 
tions in cl line-width above 30 kms -1 up to 100 kms . 
The bright background and variable extinction lends sup- 
port to the fact that an inevitable additional contribu- 
tion must result from including multiple, superimposed, 
small-scale kinematical components along every sight- 



line, some of which may be of fairly high velocity, within 
the small-scale emitting regions in the core of M82 (Sec- 
tion l5.ip . We already alluded to this effect above, where 
we invoked the unresolved nature of the line profile to 
explain the increased width of the profiles where only a 
single-Gaussian was requir ed for a satis f actory fit. In fact 
observations presented by IRozas et~a l. ( 2006) show that 
multiple emission line components are very common in 
H II regions, and can be interpreted as expanding shells 
driven by the mechanical energy input of SNe. 

A ubiquitous broad underlying component (FWHM 
~ 150-250 kms -1 ) is observed in both slits where 
the Ha line is of sufficient S/N, which cannot be 
broken down into individual subcomponents. Like 
the phenomenon of the supersonic narrow component, 
broad underlying profiles have also been observed in 
many intense star- formation sites (see in particular 
IMelnick et al] 1 19991) . but to our knowledge have never 
been studied in detail in the central regions of M82. 
Observations of the broad component in other systems 
are less well constrained, particularly for more distant 
galaxies, so the proposed explanations cited in the 
literature are more varied. These include: rotational 
effects, direct observations of stellar / cluster winds 
(|Gonzalez-Delgado et al. I1994D : the effects of SN rem- 



nants (SNRs: iRov et al.Nl992t llzotov et alj 119961): large 
scale sup erbubble expan sio n and /or blow-out (IRov et al 
1992| IMarlowe et all 119951 : iTenorio-Tagle et al 
19971 ) ; multiple resolved/unresolved expandin g shells 



(jHomeier fc Gallagher! 1 19991: IRozas et al.l l2006h: cham- 
pagne flows (lYorke et al.lll984l : lTenorio-Tagle et ai]|2000l ; 
IMelnick et al.lll999D : and ablation of gas from molecular 
clumps resulting from the impact of stellar/cluster winds 
(jMelnick et al.lfl999l ). 

Only one of these proposed explanations, that of galac- 
tic rotation, is not due to the effects of massive stars, 
and may be discounted by considering the following: the 
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maximum radial velocity relative to the systemic veloc- 
ity of M82, v™^, is of the order 100 kms" 1 . If this were 
purely caused by rotation, then, when looking through 
one side of the galaxy, the line width would also be 
~100 kms -1 . However, we see line widths of up to 
300 kms , thus proving that the broad line component 
cannot result from simple rotation. In reality, the u ra( j 
observed is due to large-scale orbital motion, not local 
gas rotation, and rotational broadening is likely close to 
negligible. 

We can also discount direct observations of the stel- 
lar winds themselves from the spatial distribution of the 
line emission. Similarly we are not directly seeing su- 
pernova remnants (SNRs) as we do not find extensive 
evidence for shock spectra. The effects of large-scale su- 
perbubblc blow-out are highly position dependent and 
would be expected to produce strong, clearly observable 
line-splitting, thus not making this a suitable explana- 
tion for the observed line broadening. Although multiple 
unresolved shells are very likely to contribute to the line 
broadening of cl, it is difficult to understand how mul- 
tiple discrete kinematical components could produce a 
smooth, Gaussian shaped, 200-300 kms -1 wide profile 
with no evidence of the discr ete components themselves 
(see also iMelnick et aT1ll999t ). Champagne flows are not 
expected to produce outflow speeds much higher than the 
sound speed in the photoionized gas (10 kms -1 for gas of 
T ~ 10 4 K), and are therefore not capable of producing 
the level of broadening we observe. Furthermore, many 
of the previous studies for which these mechanisms were 
proposed are of much larger physical regions than what 
we are observing, meaning that they may not be that ap- 
propriate to our ^10 pc scale observations. This leaves 
only the final explanation: that of hydrodynamical evap- 
oration/ablation of gas from the surface of illuminated 
gas clouds. 

These arguments are strongly supported by results 
from an analysis of spatially-resolved spectroscopy of a 
similar, high-energy environment in the starburst galaxy 
NGC 1569, presented in Westmoquette et al. (2007a, 
2007b). Here, we find striking spatial correlations be- 
tween the width of c2 and the intensity of cl, which 
strongly supports the hydrodynamic interaction mecha- 
nism. 

To explain this mechanism, we can envision a situa- 
tion whereby the energetics and dynamics of the bright- 
est clumps are driven by an interaction with the high- 
energy, ionizing photons and fast-flowing wind from the 
surrounding star clusters, and that all the observed 
emission originates from a thin interface layer at the 
clump surfaces. As the gas clouds are impacted by 
hot, fast, cluster winds, thermal evaporation and/or 
mechanical ablation of gas from their surfaces pro- 
duces a highly turbulent velocity field that pervades 
the who le region as the gas becomes entrained into 
the flow (ICharnley et alj|1990l: iBegel man fc Fabianlfl990l : 
iSlavin et all 119931 : iPittardl l2007f ) . We explore how this 
conclusion relates to the galactic wind as a whole in Sec- 
ti on 16.31 

lLord et all (|1996l ) predict that the ionizing star clus- 
ters are well mixed with molecular matter in the form 
of clouds of characteristic sizes, r c \ < 1 pc (Section [1]), 
thus providing copious sites from which gas can be evap- 
orated and/or ablated. Thus, the fact that multiple PDR 



regions are predicted from far-IR spectroscopy, together 
with the observed close packing of the individual sources, 
the uniformity of the starburst conditions throughout 
the central regions, and the identification of a ubiqui- 
tous broad- line component, all point towards gas evap- 
oration/ablation by cluster winds causing a highly tur- 
bulent velocity field being the root cause of the observed 
broad line widths. 

4. PROPERTIES OF THE IONIZED GAS 
4.1. Flux Variations 

Fig. [7] shows the distribution of the integrated Ha line- 
flux in cl and c2 along the section of the slits highlighted 
in white in Fig. [TJ Also plotted are the continuum flux 
levels (measured in the range 6620-6700 A) for each slit, 
multiplied and offset for ease of comparison (factors are 
listed in the figure caption). The relative continuum flux 
levels are consistent with the HST F814W image, and 
clearly show the bright clumps A, D and E. These re- 
gions can also be seen in both the flux levels of cl and 
c2. In region A, the two peaks in cl correspond well 
to peaks in the continuum brightness, but at the loca- 
tion of M82-A1 (offset = 0"), the Ha intensity does not 
peak so strongly. This could be due to a spatial sampling 
effect if the continuum and Ha intensities are offset by 
less than one spatial resolution element, but if the lev- 
els shown are correct, it would suggest that the ionizing 
flux emerging from the core of region A is higher than 
at M82-A1, possibly indicating that the age of region A 
is even younger than that of M82-A1 (6.5 Mvr: iPaper J . 
see also Section [2~2]) . 

The flux in c2 relative to cl indicates the significance 
of the c2 emitting gas. c2 is less than or equal to half the 
intensity of cl in the north-eastern half of region A and 
most of region C. c2 has an equal intensity to cl in the 
south-western half of region A and in region D (but here 
the peak is spatially offset by 2"). In region E, c2 peaks 
at twice the brightness of cl (but again is spatially offset 
by ^2"). A close examination of the line profiles at this 
position shows that c2 has likely been misidentified, and 
should really be assigned to c3. However, the flux peaks 
are real, and are definitely offset from the continuum 
peak. 

4.2. Interstellar Extinction 

The starburst core of M82 suffers from a high and 
non-uniform extinction. For example, several stud- 
ies have found values for the visual extinction, Ay, 
of 2-12 mag from near-IR nebul ar line ratios (assum- 
ing a foreground screen mo del) (jSatvapal et alJ 119951 : 
IForster Schreiber et al.ll2001h . In this analysis we adopt 
a foreground screen model for simplicity, although we re- 
alise that this is likely to be unrealistic considering the 
level of dust present in the M82 disc (values of Ay « 20- 
40 mag have been derived using m ixed gas and dust 
models; IForster Schreiber et al.l [20011 ). One consolation 
is that at optical wavelengths, we can only see the least 
obscured part of the starburst, thus the effect of mixed 
dust will be fairly minimal. Using slits oriented along the 
minor-axis, iHeckman et al.l (|1990t ) also measure the ex- 
tinction in the central regions of M82, and find an excess 
to the north of the nucleus up to a radius of ~300 pc, and 
that in the south, the extinction falls rapidly to a small 
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FlG. 7. — Ha flux for cl and c2 measured from the A10 spectrum sets. This time we have split the measurements for the two slits 
up, with the slit Al results in the lower panel and slit B2 results in the upper panel. The solid lines show the continuum level (in 
10~ 15 erg" 1 s" 1 cm" 2 , multiplied by 30 and offset by 800 for slit Al and by 500 for slit B2) measured in the range 6620-6700 A. (A color 
version of this figure is available on-line.) 

and constant value beyond R « 300 pc. These findings 
are consistent with the assumed inclination of the galaxy. 

To calculate the reddening distribution of the gas we 
compared the observed flux ratio of H/3/Ha (assuming 
no underlying stel lar absorption) to t h e the oretical case 
B ratio given by iHummer fc Storeyl (11987D , using the 

th| (|1983l ). Shown in 



Galactic extinction law of iHowartl 
Fig. [5] is the extinction in magnitudes measured at 20 
pixel (1") intervals along the slits. The uncertainty in 
the measurement of the H/3 line flux dominates the er- 
ror bars shown, and indeed only four extracted spectra 
from slit B2 had detectable H/3 emission. For the few 
extracted spectra from slit Al where we detect H7, we 
also derived Ay from the Hy/H/3 ratio. The results were 
consistent with the extinction values obtained using the 
Ha/H/3 ratio, although the errors were large. 

The average extinction in region A is Ay = 5.5 mag, 
which falls to an average of 4 mag in region C, and the 
data are consistent with a further decrease past region C 



towards the south-west. The extinction towards M82-A1 
(|Paper II ) is less than that towards the majority of re- 
gion A, perhaps explaining why this cluster stands out so 
clearly in optical images. We can only measure the red- 
dening for one point to the north-east of M82-A1, and 
although the errors are large, it is consistent with be- 
ing at a much higher value than the surroundings. This 
coincides with the beginning of a region of complete ob- 
scuration at optical wavelengths, extending towards the 
north-east of the starburst core (the dust-lane) . Extinc- 
tion measurements of the diffuse gas in this region (—10 
to —2 arcsecs) are therefore made impossible due to the 
very low S/N or absence of the H/3 lines. 

The extinction in slit B2 could only be measured for 
four points representing the brightest emission in the 
cores of regions D and E. In region D we find Ay « 
4 mag, and is lower than that in region A. In region E, 
however, the extinction appears to be much higher at 
Ay w 6-6.5 mag. 
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4.3. Electron Density 

We have measured the [Sn]A6717/A6731 line ratio by 
fitting Gaussian profiles to the emission lines in each of 
the A20 spectra for both slits. The distribution of de- 
rived electron densities, n e , (assuming an electron tem- 
perature, T c = 10 4 K) is plotted in Fig. [9] 
In slit Al, we find the electron densit y to be highest 



around M82-A1 



isoolSjg 



*: IPaper J ). [Sn] 

emission is not detected at all in spectra extracted from 
the north-east of M82-A1 (negative offsets) , so in an at- 
tempt to derive an approximate average density for this 
region, we extracted a 60 pixel wide spectrum. How- 
ever even summing over this many pixels did not yield a 
[S n] detection, and therefore no density could be derived. 
From M82-A1, the density falls rapidly to ~1000 cm" 3 in 
the centre of region A and to 600-800 cm -3 in region C. 
Our observations provide constraints on the steepness of 
the density gradient from region A to Al: an enhance- 
ment of n c ~ 1000 cm -3 takes place over a projected 
distance of only 2.5 arcsecs (45 pc). In the area located 
between the two cluster complexes A and C, n c falls be- 
low the low-density limit (<100 cm -3 ) for the [Sn] ratio 
method, and is indicated by upper-limit arrows in Fig. [9j 
Our measurements are lower than those of 
lO'Connell fc Manganol (|1978f ). who also measured 
n e from the [Sn] doublet ratio for regions A and C, and 
found an average density of n Q — 1800 cm" 3 . The dis- 
crepancy may arise because at lower spatial-resolution, 
their data inevitably suffered from a higher degree 
of contamination and luminosity weighting than our 
high-resolution data, and were therefore biassed towards 



the densest gas. iHeckman et al.1 (|1990| ) also measured 
the electron density from their minor-axis slits, and 
found an n c of ~1000 cm -3 at the nucleus, falling to 
~100-300 cm" 3 at lateral distances of 1-2 kpc, which is 
in good agreement with our measurements. 

We measure a very different distribution of densities 
in slit B2: close to the centre of region D 8 , the electron 
density is as low as 300 cm -3 , but quickly rises to over 
2000 cm" 3 in the area between regions D and E. Here 
slit B2 passes just to the north of a smaller cluster com- 
plex (2-4 arcsecs north of the 2.2 ^m nucleus; Fig. [1]), 
which if really surrounded by gas at this density would be 
a very stifled region indeed. Alternatively, this density 
peak could result from scattered light from the starburst 
core/nucleus, under the assumption that higher pressures 
are associated with more intense zones of star formation. 
In region E, the electron density is also high with val- 
ues of ^1300-1700 cm" 3 , similar to that measured near 
M82-A1. 



5. PHYSICAL STATE OF THE IONIZED ISM 

5.1. Emission Measure 

In IPaper ll we discussed the surprising uniformity 
of the starburst over a large scale range by compar- 
ing our calculation of the ionization parameter of the 
M82-A1 Hn region to the r esults of the IR study by 
iForster Schreiber et al.l (|200lh . We now use the observed 
Ha fluxes to derive an estimate of the scale-length of the 

8 Note: slit B2 is not aligned with the major axis of region D 
and misses most of the eastern extent. 
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Fig. 9. — Electron density derived from the [S II] A6717/A 6731 flux ratio, measured from the A20 spectrum sets. The triangle represents 
the electron density of the Hn region surrounding M82-A1 IjPaper HI . and the arrows represent slit Al measurements which fell below the 
low density limit for the [S II] method. 

emitting regions to compare to and extend this discus- 
sion. 

A characteristic scale-length can be derived from the 
emission-measure, which is defined by: 



pint 



2 Hq hvQlot) 
ni ogg A _^, 9 e dV, 



AttD 2 



(1) 



where a^g is the case B recombination coefficient for 
Ha (8.64 x 1CT 14 cm 3 s" 1 : IOsterbroclH98l . hv(Ra) is 
the energy of an Ha photon, D is the distance to the 
source, and e is the gas volume filling factor. To remove 
the distance dependence, and to re-arrange in terms of 
a scale-length, dV can be substituted with dA d£, where 
dA = A/AttD 2 = angular area on the sky (0.1 x 1 arcsec 
in this case). Integrating over the whole emitting volume 
gives the mean intensity of Ha, 



erg s 1 cm 2 arcsec . 



= K>«& (e)£ 

Hence, in reality what we can measure is £(e = 1) for 
all three Ha line components in slit Al. Since e < 1.0, 
this is a lower-limit to the size-scales of the emitting H II 
columns. 

We find £(e = 1) to be less than 100 pc at all points, 
and is as low as ~few parsecs to a few tens of parsecs in 
the north-east half of region A and the majority of re- 
gion C. The consistency between line components shows 
that the gas from which each is emitted has similar char- 
acteristic sizes. Within the south-western half of region 
A and the A-C inter-clump region, the scale-lengths are 
larger, rising up to ~50 pc. 



The roughly constant values of the ionization parame- 
ter, log U « —2.3, within the main M82 starburst zone 
suggests that dust could be successfully competing with 
gas for ionizing ph otons on the sur faces of cloud s (see 
iDopita et~ al. 2002). In this case the lDopita et al.l model 
also suggests radiation pressure could be a significant fac- 
tor. The luminosity densities in the M82 starburst zone, 
however, are well below those of the narrow line regions 
of AGN (e.g.lCecil et al.l l200l) or found in ULIRGs (e.g. 
iThompson et al.l I2005Q where models show that radia- 
tion acting on dust becomes important. Consistent with 
this view, we also find that the energy density in turbu- 
lent ISM motions exceeds the thermal pressure, and so 
even if radiation and thermal pressure are comparable, 
they are unlikely to be strongly influencing the M82 gas 
dynamics. 

5.2. Comparison to ISM models 

We can now discuss how our findings relate to the ISM 
models introduced in Section [TJ and what that implies 
for the M82 starburst. In the above section we found 
that the characteristic scale-size of the emitting regions 
remains below ~50 pc (3") throughout regions A and C. 
The uniformity of the starburst conditions must there- 
fore stem from the fact that the scale of t he emitting 
regio ns remains very small (as suggested by lLord et all 
Il996| ) , hence the conditions are able to respond rapidly in 
order for the starburst to maintain such a near-constant 
state. However, within this overall uniformity there are 
variations. The smallest characteristic sizes (~few par- 
secs) are found in or near the dense cores of the clumps, 
whereas the largest sizes arc found in the inter-clump re- 
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gion. These differences imply that the most compact 
clouds are found in the clump cores where the star- 
formation is most intense and the gas pressures are high. 

Diagnostic flux ratios derived from the most com- 
monly observed nebular lines can give a further indi- 
cation of the ionization parameter and the mean level 
of ionization within the M82 starburst. The ratios 
of [Sn]A6717+A6731/Ha and [N n]A6548+A6583/Ha are 
primarily tracers of the ion ization parameter of the gas 
(jDopita et ai]|2000t I2006H) . and are shown in Fig. [TO] 
as a function of position along slit Al. We find that 
the [S ii] /Ha ratio changes very slowly across regions A 
and C, from ~-1.8 dex at M82-A1, to ~-l dex to the 
south-west of region C. For [Nii]/Ha, we were able to 
measure the ratio independently for both cl and c2 in 
many places, but find no significant difference between 
the two (except for perhaps in the inter-clump region 
where c2 is consistently higher than cl for seven con- 
secutive points). Again, the ratio stays roughly constant 
across regions A and C at an average of — 0.2±0.5 d ex. As 
a comparison, iShopbell fc Bland-HawthoTnl (fl998h mea- 
sure equally high [Nii]/Ha ratios (—0.5 to —0.2 dex) in 
the inner collimated zone for their low velocity compo- 
nent. Having access to the bluer nebular lines of [Oin] 
and H/3 means that we can use their ratio as an indica- 
tor of the mean temperature (r e ff) of the ionizin g sources 
(|Veilleux fc Osterbrockll987l : lDopita et al.ll2000l ). Unfor- 
tunately though, our ability to measure this ratio, like 
that of measuring the extinction, is limited by the low 
S/N of the H/? line. Where this line is detected, we find 
the flux ratio to remain constant within the uncertainties 
at log([Oili1/H/3) pa —0 .5 over both regions A and C. 

IRigbv fc Riekd (|2004f ) discuss the effectiveness of cer- 
tain optical and IR nebular line ratios, including 
[O in] /H/3, at tracing the mean temperature of the ion- 
ization source by comparing the consistency of the indi- 
cators between themselves and to models. They suggest 
that finding r atios that imply a low mean T B g in starburst 
galaxies (e.g. iThornlev et al1l2000l . this work) could re- 
sult from either a lack of massive stars (M > 40 M Q ), 
or more likely, from the fact that the ultracompact Hn 
region phase is prolonged by the high densities and pres- 
sures in starburst ISMs. Under the influence of the high 
pressures, the most massive stars within the clusters may 
spend a significant proportion of their lifetimes embed- 
ded within the den se, highly-extincted (Ay < 50 mag; 
lHanson et aT1l2002[ ) regions in which they form, and they 
would remain undetectable, even through re-processed 
nebular optical or IR emission. Although their bolo- 
metric output may contribute to the surrounding gas for 
their entire main-sequence lifetime, it may influence neb- 
ular excitation indi cators for a smaller fraction of this 
time (~85 per cent; IRigbv fc Riekdl200l . 

This situation agrees with what we know about the 
environment in M82, including its remarkably homo- 
geneo us ionization parameter (|Forster Schreiber et al.l 
f200lh and our detectio n of the dense, compact Hn region 
surrounding M82-A1 (jPaper J ). Clearly not all clusters 
can be deeply embedded, since in many cases we can 
actually see the starlight from the cluster directly. How- 
ever, the consistency in the line ratios (tracing nebular 
excitation and T c g) across the starburst regions must re- 
sult from both the highly fractious nature of the ISM and 
the uniformity of the escaping radiation field. This radi- 



ation field must be attenuated by the high density gas to 
an extent to which it no longer reflects the presumably 
very inhomogeneous distribution of high-energy photons 
from the most massive stars within the clusters. 

6. STRUCTURE OF THE M82 STARBURST 

6.1. Evidence of the bar 

We have used the high spectral- and spatial-resolution 
of our observations to track the kinematics of the Ha line 
components across the two slits in unprecedented detail. 
We now use these measurements to provide constraints 
on the orbital parameters of the ionized gas, and the 
location of the individual starburst clumps within the 
galaxy. 

As discussed in Section [Q IWills et all <|2000D devel- 
oped a model for the M82 bar system by analysing ve- 
locities derived from a number of neutral and ionized 
emission lines. They found the best-fitting bar model 
to have a peak radial velocity projected onto the sky of 
140 kms -1 , an angular velocity of 217 kms" 1 arcsec -1 , 
a total length of 1 kpc, a core radius of 25 pc, but as- 
sume an opposite inclination to that which is commonly 
accepted (i.e. —80°). We now compare the results of 
their model to our observations to determine if we sec 
evidence for any of the orbit families. In Fig. QTJ we 
plot the observed major-axis Call stellar absorption and 
[Sih] A9069 and Pa( 1 0) ne bular emission-line velocities 
from iMcKeith et alJ l)1993f) and ou r Ha cl and c3 ve- 
locities (see also iGreve et al.l l2002f ). The a^-orbits, as 
traced by the stellar Can measurements, can clearly be 
seen departing from the velocities determined from the 
gas emission-lines at offsets of ^+13" and ~0". Interest- 
ingly, M82-A1 appears to be located near the very top of 
the X2-orbit pattern, unambiguously separated from the 
Xi-orbits. In other parts of the plot, the close match be- 
tween the IR and optical data-sets indicates that the Ha- 
emitting gas follows the stellar rotation closely, although 
the rapidly varying extinction introduces inevitable com- 
plications. The inflection points at ~0 and ~+17" repre- 
sent the intersection between the x<i- and a?i-like orbits, 
where the gas is expected to shock and begin its journey 
towards the nucleus. Region C is clearly not associated 
with the X2-like orbits, implying that it may be located 
within the molecular ring of material beyond the bar. 

6.2. Locating the starburst zones 

Bearing in mind that the IWills et al.l bar model as- 
sumed the opposite inclination to that which is gener- 
ally accepted, we can now re-asses how the bar might 
be oriented with respect to the line -of-sight, and to the 
individual clumps. The lLord et alJ ISM model predicts 
a concentration of molecular clouds (PDR regions) at ei- 
ther end of the stellar bar, which may be attributed to the 
interaction of bar x\- and a^-orbits causing the build-up 
of dust along the leadin g xi-orbits (a phenomenon some- 
times termed 'spraying '• |Athanassoulall992rirWills et all 
12000( 1 . Fig.[12lsho ws contours of inte grated CO (J=l->0) 
line intensity from lWeifi et alJ (|2001( 1 overlaid on an HST 
broad-band colour-composite, clearly showing how the 
peak of the CO emission corresponds to the location of 
the dust clouds either side of the nucleus a nd the pro- 
posed position of the two PDR hot-spots (|Lord et al.l 
11996ft . 
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Fig. 10.— Plots showing the variation in the flux ratio of [S n] A6717+A6731/Ha (left) and [N n] A6548+A6583/Ha (right) with position 
on the Al slit. The S/N of the Ha and [Nil] lines are high enough for the excitation of the two line components to be shown separately. 
For comparison, the continuum flux level is plotted as a solid line (arbitrary scale) on both plots. 
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Fig. 11. — Ha cl and c3 radial velocities for slit Al plotted 
together with the near-IR stellar absorption line Ca II A8542 and 
the [S III] A9069 and Pa(10)A9014 emission lines from lMcKeith et al.l 
(119931 ') . Offsets are measured relative to the position of M82-A1. 




Fi g. 12. — Contou rs of integrated C ls O(J=l— >0) line inten- 
sity IWeifi et al.l 120011) overlaid on an HST broad-band colour- 
composite (blue: F439W, green: F555W, red: F814W; credit 
NASA/ESA R. de Grijs). The CO peaks trace the molecular /dust 
torus surrounding the stellar bar on either side of the nucleus. The 
location of the two STIS slits are marked with solid lines; the ma- 
jor axis is indicated with a dashed line and the cross marks the 
position of the 2.2 lira nucleus. 

We can therefore associate this pile-up of molecular 
material/dust at either ends of the bar with the highly 
obscured regions to the north-east of clumps A and D 
and between clumps C and E extending out to the south- 
west. Our me asurement of a rapid increase in extinction 
(also found bv lSatvapal et al.lll995HAlonso-Herrero et al.1 
I2003T ) and electron density to the north-east of regions A 
and D (see also Fig. [T2")) . and the high extinction and den- 
sity in region E and the D-E inter-clump region, is con- 



sistent with this conclusion. Evidence for a rapid fall-off 
in extinction to the south-west of region C may indicate 
the torus edge in this direction (Fig. [5]). The size and 
location of these dust clouds is consistent with the pre- 
dicted bar length (~1 kpc), and the comparative levels 
of extinction and prominence of the north-eastern dust 
cloud (again see Fig. [T2")) would imply that this is on the 
near-side of the starburst. This dust cloud also appears 
to wrap around region D and extend up into the north- 
ern inner- wind region - whether or not this dust is being 
entrained into the wind flow at this point would require 
detailed kincmatical observations. 

In the region between clumps A and C, we do not mea- 
sure a discontinuity in the extinction, but we do find a 
sudden drop in the electron density (Fig. [9|) . This sug- 
gests that there is a physical 'hole' in the gas distribution 
between these two clumps and that the two regions are 
physically distinct, rather than the decrease in surface- 
brightness being due to an obscuring dust-lane. The re- 
maining gas in this 'void' must still be kincmatically as- 
so ciated with the two clumps ( confirming the suggestion 
of lO'Connell fc Manganol[l97^ since there is no velocity 
discontinuity here. We can therefore now associate re- 
gion C, located to the south of the bar region, as being 
outside the X2-orbit region, and part of (or a result of 
star-formation within) the molecular torus surrounding 
the bar. 

The fact that we can see evidence for both ends of the 
bar indicates that we are able to sample light from a con- 
siderable fraction of the M82 system. Contrary to previ- 
ous models of the M82 system that assume la rge fractions 
of the M82 system are optic ally invisible (e.g. lRieke et al.l 
119931 ISatvapal et al]fl~997t ). the main hidden region ap- 
pears to be only in the middle of the bar near the nucleus. 
The most recent and perhaps most intense star-forming 
sites are located inside compact clouds spread through- 
out the system but concentrated in the torus surrounding 
the bar. 

The fact that M82-A1 appears to be located at one end 
of the X2-orbits places a strong constraint on the orien- 
tation of the bar within M82, and indicates that region 
A may have formed as a result of intense star-formation 
in the intersection between the x\- and X2-orbit families. 
The shallower velocity gradient between regions D and E 
indicate that they are located at a larger radius than that 
of regions A and C, but the distinction between these two 
regions is less clear than that between clumps A and C. 
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Fig. 13. — Our proposed orientation of the M82 bar and disc, with respect to the main starburst clumps (labelled A, C, D and E) within 
the inner 2 kpc of the galaxy. The disc is inclined 80° to the line-of-sight and at a PA of 70°. The bar xi- and X2-orbits are shown with 
an outlined and filled ellipse, respectively, and are also reproduced in plan-view in the upper-left. The major-axis of the bar extends along 
the :ri-orbits and is ~1 kpc in length. At the ends of the bar, molecular material and dust (hatched regions) correspond to the location 
of the PDR hot-spots in the ILord et al.l J1996I ) ISM model. The cluster M82-A1 (shown as a white plus-sign) is located at the end of an 
X2-orbit with the maximum redshift observed. Clump C shows the highest blucshift indicating that it is located at the opposite end of 
the bar, whereas regions D and E must be located either behind or in front of the nucleus (shown as a white cross) due to the shallower 
velocity gradient observed be twe en the two. The wind outflow cones are indicated above and below the disc. For comparison, it may be 
helpful to refer to Figs [T] and IT2l 



Comparing Fig. Q] to Fig. [12] shows that there is a great 
deal of extended emission and variable extinction in this 
inter-clump region, and the confused and variable veloci- 
ties in this region indicate that we are seeing glimpses of 
star clusters through holes in a foreground dust screen. 
However, the fact that we find narrower line widths in 
regions D and E compared to A and C (Section l3?2|) sup- 
ports the scenario of them being located at larger radii, 
since narrower line widths would indicate a more qui- 
escent gas state, consistent with material further away 
from the nucleus. 

Fig. [121 shows a schematic of our proposed spatial con- 
figuration of the starburst clumps and dust lanes in re- 
lation to the orientation of the bar orbits. This model 
implies that the eastern side of the bar is more distant, 
and that it extends from behind the dense dust cloud 
to the east of region A to just past region E. This ori- 
entation also implies that the eastern side of the inner 
X2-orbits, where M82-A1 is located, is nearest. 

6.3. Implications for the production of the galactic wind 

In this section, we will use our observations together 
with the discussion given in the previous two sections to 
consider what implications our findings may have on how 



the superwind is produced in the starburst core. We will 
first consider how our spatially-resolved measurements of 
the gas density may affect the individual cluster winds, 
then go on to discuss possible explanations for the width 
of the observed line profiles, and finally how the dynamics 
of the broad component may indicate that we are seeing 
the roots of the large-scale galactic superwind. 

6.3.1. Densities and pressures 

In lPaper"H , we explored how the unusually high density 
of the gas surrounding M82-A1 has stifled the growth of 
its Hii region, and commented how it is hard to under- 
stand the development of the superwind if this cluster is 
representative of the entire starburst. We are now in a 
better position to know what the representative condi- 
tions are, and how these might vary over the face of the 
starburst. It is clear from high- resolution HST images 
of the galactic wind (Gallagher et al., in prep.) that the 
outflow is being driven directly from the central clumps 
(region A in particular), but our measured densities of 
the gas within these cluster complexes arc still very high 
(500-900 cm" 3 , or P/k « 0.5-1.0xl0 7 cm" 3 K; lPaper J ) 
compared to standard Galactic ISM values. However, 
they are less than half of that found near M82-A1, and 
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this highlights a possibly significant difference between 
a relatively isolated cluster such as M82-A1 and the 
densely-packed clusters in the clump cores. It appears 
that instead of having collimated winds from individ- 
ual SSCs, the cluster complexes (clumps) produce a high 
pressure zone, or an 'energy injection centre', and the 
winds result from hot gas expanding out of these zones, 
upwards into the halo. This scenario will be explored in 
more detail by Gallagher et al. (in prep.). 

6.3.2. Dynamics of the broad line component 

We can now turn our attention to whether, by observ- 
ing this broad, underlying component to the Ha line, 
we are indirectly detecting the roots of the supcrwind 
flow. In our proposed hypothesis, c2 represents gas that 
has been stripped off the surface of gas clouds embed- 
ded in the starburst, so in principle the radial velocity 
offset between this component and the brighter (hence 
presumably denser) cl gas can be used to trace the flow 
of the hot superwind into which it is being entrained. 

Slit Al passes through the very core of region A, 
and even though here we see some of the broadest line 
widths (up to 300 kms" 1 ), we do not see any velocity 
offset between the components, indicating that at this 
point there is still no preferred outflow direction. In 
region C, the slit passes through the southern side of 
this clump, rather than its core. Already at this ra- 
dius, some evidence of an or dered wind flow is det ected 
in the kinematical studies of IMcKeith et "aTl (|1995| ) and 
IShopbell fc Bland-Hawthornl (|1998fl . The"se~i authors find 
that within 200 pc south of the nucleus the Ha com- 
ponents are separated by ~50 kms" 1 . However, this 
is still within the chaotic zone of the outflow, and cer- 
tainly well within the energy injection zone, where it is 
thought th e wind flow is still subsonic but r apidly ac- 
celerating (jShopbell fc Bland-Hawthorn! [19981 ). In this 
region we find that c2 is redshifted with respect to cl 
by ~40 kins" 1 and we also detect an ~80 kins" 1 red- 
shifted c3 profile (see Fig. [5|. These values are consistent 
with the aforementioned studies, and show that here we 
are beginning to see an ordered flow with a preferred 
direction, and that c2 is tracing the wind roots. 

The region beyond clump C is interesting in a num- 
ber of ways: here the cl velocities rapidly decrease back 
to w sys , d eparting from t h e maj or-axis velocity measure- 
ments of IMcKeith et all (|1993| i. and the offset between 
cl and c2 reverses direction, becoming blucshifted by 
~40 kins" 1 relative to cl. This far end of the slit is co- 
incident with lateral ionized streamers pointing towards 
the south-west (see Fig. Q]), indicating that this region 
may be part of a distinct flow, separate to the main 
wind. If clump C is located within the molecular torus 
surrounding the bar, then this region may be less in- 
fluenced by the bar's collimating effects, explaining the 
lateral flow. 

7. SUMMARY 

The use of high spatial resolution HST optical images 
and STIS spectra have enabled us to fill in a few more 
pieces of the large and complex jigsaw puzzle of M82. 
The data-set has allowed us to probe the ionized gas en- 
vironment across the nuclear starburst in unprecedented 
detail. We have concentrated on mapping the ionized 
gas properties in the starburst core using two STIS slit 



pointings. One slit passes through the core of region A, 
and the southern edge of region C, whilst the other slit 
passes through part of region D and ends within region 
E. 

• By performing a comparison of our radial velocity 
measurements to previous observations and mod- 
els, we confirm the presence of a stellar bar with 
an inner Lindblad resonance (ILR), resulting in a 
set of unique family of .T2-orbits perpendicular to 
the bar major-axis. The radial velocity of M82- 
Al is consistent with being on an X2-orbit at one 
end of the perpendicular ellipse, implying that re- 
gion A may have formed as a result of intense 
star-formation at the intersection o f the x\- and 
X2-orb its (as predicted by models; lAthanassoulal 
Il992ah . 

• Since the interaction of xi-orbits with the perpen- 
dicular X2-orbits is predicted to lead to a build- 
up of gas and dust along th e leading x\ -orbits 
(lAthanassoula fc Bureau! fl999T ) . we conclude that 
the prominent dust clouds to the north-east of re- 
gions A and D, and between regions C and E ex- 
tending towards the south-west, must have been 
formed in this way. 

• By comparing CO observations to an HST broad- 
band colour composite of M82, we identify that 
the peak of the molecular emission corresponds 
to t he molecular torus surrounding the ionize d 
ring (jAchtermann fc Lacvlfl995t IWeifl et al.ll200lh . 
the propos e d pos itions of the PDR hotspots from 
lLord et all (|1996l ). and the prominent dust clouds 
located on either side of the nucleus. Since region 
C is clearly not associated with the a^-orbits, we 
assume it to be located within this molecular torus. 

• Our data are consistent with a physical distinc- 
tion between clumps A and C, but gas between 
the clumps follows the same radial velocity pattern 
as the clumps, suggesting that all of this material 
is in the rotating disk of M82. 

• A shallower velocity gradient between regions D 
and E leads us to believe that they are located 
at a larger radius to clumps A and C, supported 
by the finding of narrower line widths in these re- 
gions. However, the strong, patchy foreground ex- 
tinction indicates that we may be seeing these clus- 
ter complexes through gaps in a thick foreground 
dust screen. 

• By measuring the [S n] doublet ratio along the slits, 
we find the gas densities (pressures) in the individ- 
ual clumps to be high (500-900 cm" 3 ; P/k 0.5- 
1.0 x 10 7 cm -3 K), but significantly lower than 
that measured for the gas immediately surrounding 
M82-A1. We find that instead of seeing collimated 
winds from individual SSCs, the cluster complexes 
(clumps) each appear to produce a high pressure 
zone, or an 'energy injection centre', and that the 
winds result from hot gas expanding out of these 
zones. We will explore this scenario in more detail 
in a forthcoming contribution (Gallagher et al. in 
prep.). 
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• We exten d the uniformity of the ISM conditions 
found by iForster Schreiber et al.1 (|200lh down to 
parsec-scales by deriving the emitting scale-sizes 
and by measuring the ionization parameter sensi- 
tive line-ratios of [Sn]/Ha and [Nil] /Ha. We find 
the size-scale of the emitting regions to be always 
less than 50 pc along slit Al, and that the small- 
est regions (few pc) are in or near the dense cores 
of clumps A and C, implying that the most com- 
pact clouds are found where the star-formation is 
most intense. Our results confirm the prediction 
of uniform conditions throughout the starburst re- 
sulting from the highly fragmented nature of the 
ISM and/or the importance of dust in absorbing a 
significant fraction of the Lyman continuum lumi- 
nosity. Our data also support the well-mixed na- 
ture of the embedded ionizing sources and ionized 
gas. 

• Our line-profile analysis has shown that the widths 
of all but a few line components are >30 kms -1 . 
We have explored the possible mechanisms for pro- 
ducing the observed broad-line widths in both the 
narrow component (cl) and the broad component 
(c2). The most likely explanation for broadening 
of cl is a combination of gravitationally induced 
virial effects and the stirring of the ISM through 
the intense star-formation activity. A base level 
of turbulent broadening (^30 kms™ 1 ) exists over 
the whole starburst region, whilst the scatter to 
broader widths (up to 100 kms" 1 ) is likely to re- 
sult from the presence of multiple unresolved kine- 
matical components (e.g. expanding shells) along 
the line-of-sight. The contribution from unresolved 
components is inevitable due to the presence of 
large numbers of densely packed SSCs. 

• A ubiquitous broad component (150-250 kms -1 ) 
is seen in all nebular lines with a high enough S/N. 
We argue that the broad line widths are produced 
by the interaction of high-energy, ionizing photons 
and fast-flowing winds from the star clusters with 
the cool gas clumps found throughout the star- 
burst zone. As the wind impacts the surface of the 
clouds, evaporation and ablation of cloud material 
results in a highly turbulent velocity field which 
manifests as this broad component. 

• In our model, the broad component represents tur- 
bulent gas stripped from the surface of gas clouds 
through the action of stellar winds and SN from 
the surrounding SSCs, so therefore could be used 
to track the hot superwind as the cooler material 
becomes entrained into the wind flow. We find no 
velocity offset between cl and c2 in of region A, 
implying that here in the clump cores, the wind has 
yet to possess a preferred outflow direction. Where 
the slit passes through the southern edge of region 
C, we begin to see an ordered flow: c2 has a con- 
sistent rcdshiftcd offset, and a further rcdshiftcd 
narrow component (c3) can be identified. 



• Beyond clump C, the offset between the cl and 
c2 velocities is reversed, and cl departs from the 
radial velocity pattern measured along the major- 
axis. Since these measurements are coincident with 
a number of lateral south-west pointing streamers, 
this has lead us to believe that this region may 
not be associated with the main collimated flow. 
If clump C is located within the molecular torus 
surrounding the bar, then this region may be less 
influenced by the bar's collimating effects, explain- 
ing the lateral flow. 

In conclusion, it seems a number of key ingredients are 
needed for the M82 starburst to exist. A central concen- 
tration of gas is required, together with clumps of dense 
star formation and high interstellar pressures. Of course, 
these three essential components do not exist in isolation, 
they link together and are to a larger or lesser extent the 
result of external influences. Large-scale galaxy encoun- 
ters (which we know M82 has experienced) very often 
result in the formation of a bar, and a bar is a very ef- 
ficient method to funnel gas into the nuclear regions in 
order to build up the required central gas concentration. 
Furthermore, both clumps and bars can be the products 
of disc instabilities that occur when discs are too dynam- 
ically cool i.e. gas rich. We can therefore reason that the 
properties of the M82 starburst are all inextricably linked 
to one another. The study of one will always require a 
study (or at the least a good understanding) of the oth- 
ers. 

Furthermore, we can now explain both the highly frag- 
mented nature of the ISM and the near-uniform star- 
burst conditions through the high densities and pres- 
sures within the starburst core, since high ambient den- 
sities result in compact, small, high density interstellar 
gas clouds. In high-pressured galaxies such as M82, the 
characteristic scale-sizes are therefore smaller than less- 
pressured systems, and this means that the gas condi- 
tions are able to respond quickly to any changes. This 
ability to respond rapidly results in the observed near- 
constant state of the starburst (seen though indicators 
such as ionization parameter and line diagnostics). 
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